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ABSTRACT 

A study was conducted to evaluate the accuracy of a commercial computational fluid 
dynamics (CFD) code (CFDRC-ACE+) for predicting incompressible air jet flows with 
simple geometries. Specifically, the axis-symmetric and two-dimensional heated air-jets were 
simulated using a standard k-e turbulence model. These CFD predictions were directly 
compared to an extensive compilation of experimental data from archive literature. The round 
jet results indicated that the code over-predicted the velocity-spreading rate by 24% and the 
temperature-spreading rate by 29%. In addition, the centerline velocity and temperature decay 
rates were also over-predicted by 21% and 30%, respectively. The geometric and kinematic 
virtual origins were over-predicted, as well, by approximately 7.5 diameters for the velocity 
profiles and 10.5 diameters for the temperature profiles. The planar jet simulation was 
generally closer to experimental data ranges, with an under-prediction of the velocity- 
spreading rate of approximately 17% with an over-predicted temperature-spreading rate of 
12%. The centerline velocity and temperature decay rates were both under-predicted at 22% 
and 27%, respectively. Again, the geometric and kinematic virtual origins were over- 
predicted by approximately 7.5 slot heights for the velocity profiles and 10.5 slot heights for 


the temperature profiles. 
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I. INTRODUCTION 


A. MOTIVATION 


The characteristics of jet development are very important in many military and industrial 
applications including eductors, thrust augmenting ejectors and exhaust suppression units. For 
example, current efforts are underway to improve present shipboard exhaust eductors without 
adversely effecting other parts of the ship. Former redesign of these types of systems has mainly 
been accomplished through extensive scale model testing of each new design proposal. The 
desire to reduce developmental costs of design test configurations has stimulated an increased 
interest in Computational Fluid Dynamics (CFD) design studies. By using CFD analysis, only 


the most favorable design modifications warrant further experimental testing. 


B. CHALLENGE 


CFD simulations help researchers to understand and visualize the complex flow patterns 
before costly experiments are performed. However, CFD codes are typically challenged when 
predicting flows that are dominated by turbulent shear mixing. The general conservation 
equations of continuity, momentum and energy can be solved directly by using Direct Numerical 
Simulation (DNS) on supercomputers. However, the extremely small grid-cell size and time 
steps required to resolve the smallest scales of turbulence (Kolmogorov scale) limit DNS 


calculations to low/moderate Reynolds number (Re) flows only. 


In order to reduce the amount of computational resources needed, turbulence models 
have been developed to mimic the effects of turbulence (Reynolds stresses) that are contained 
within the conservation equations. This approach permits the use of smaller computers to model 
complex and high-speed flow configurations. However, turbulence models merely approximate 
the Reynolds stresses and are usually limited to a small number of simple flow regimes. Often 
‘turbulence model constants’ must be adjusted to make the model accurately predict other types 


of flows. 


Depending on the turbulence model used, any number of empirically derived constants 
may be required to properly model the flow. Overall, the process of correctly modeling the 
Reynolds stresses of a flow configuration is extremely difficult and problematic. Ideally, a 
calculation derived from the flow geometry itself would be the most desirable option when 
calibrating a turbulence model. However, only empirically derived constants from relatively 
simple experiments have been able to provide such inputs. Unfortunately, the highly complex 
flow patterns that typically arise in nature have precluded these experimental constants from 
being truly universal. Therefore, a “trial and error” methodology is often employed to adjust the 


turbulence model constants for reasonably accurate predictions in complex flows. 


Commercial codes, designed to model many types of flows with average operator skill, 
often make certain concessions to provide acceptable overall predictions. Frequently, the 
constants are an integral part of a commercial code and can not be changed to fit the flow regime 


under consideration. On the other hand, ‘wall functions’ and other additional turbulence features 


are commonly used to modify or correct the local flow conditions as necessary. Overall, the 
process of using turbulence models is highly challenging and the results are sometimes less than 


desired. 


C. OBJECTIVE 


To eventually predict complex multiple-jet type flows, a single free jet must first be 
correctly modeled without the interference of walls or adjacent jets. Accurate simulation of the 
simple ‘free jet’ is crucial for the continued use of the CFD code and earblenoe model. With 
this baseline of code experience and knowledge, further research can then progress towards 


developing and configuring more complex flow configurations. 


The objective of this research is to determine the performance and sensitivity of CFD 
Research Corporation’s numerical analysis code (CFD-ACE+) to predict simple nearly 
incompressible free turbulent jet flows. In particular, heated air axisymmetric (round) and two- 
dimensional (planar) jets were modeled and compared to tabulated experimental data. In 
addition, a critical review of experimental data available in the archive literature is also provided. 
The difficulty in predicting these simple flows has proven to be one of the more challenging 
problems faced by many turbulence modelers over the last five decades. The CFD-ACE+ code 
experience and jet behavior knowledge provided by this research was meant to create a solid 


foundation for further investigations. 
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II. BACKGROUND ON FREE JETS 


A. PREFACE 


The turbulent jet has been the subject of many extensive studies, both theoretical (i.e., 
Abramovich 1963, Hinze 1975, Townsend 1976, Schlichting 1979, Schetz 1980, List / Rodi 
1982), and experimental (i.e., Hinze / van der Hegge Zijnen 1949, Wygnanski / Fiedler 1969, 
Rodi 1975, Panchapakesan / Lumley 1993, Hussein et al 1994). Being both basic and important 
in many applications, ‘simple geometry’ jets are actually quite complex and multifaceted. From 
the earliest investigations, the turbulent jet has been dissected into different discrete regions for 
in-depth analysis. A brief description will be given here to familiarize the reader with this type 


of flow. 


A free jet is formed when a fluid is discharged from an orifice or nozzle into quiescent, 
co-flowing, or counter-flowing environments. Co-flows and counter-flows are not considered in 
this investigation. By analyzing only subsonic flows with Mach Numbers less than 0.3, 
compressibility effects can be ignored. In addition, only jets with “top-hat” velocity / 


temperature exit profiles, negligible buoyancy effects, and stagnant surroundings are considered. 
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Figure 1. A Typical Free Jet 


As shown in Figure 1, the jet is usually divided into three main regions. Each region has 
its own unique characteristics that set it apart from the others. The virtual origin is 
conceptualized as an “ideal” jet source that produces the same self-similar jet effects without the 
transient regions (potential core and intermediate). The “geometric” virtual origin, as indicated 
in Figure 1, is obtained by extending the linear spreading profile of the jet backward to the jet’s 
centerline. The distance from the jet center (where the velocity is Uc) to the half-velocity point 
(Uc/2) usually defines the jet’s spreading profile length. The “kinematic” virtual origin, on the 
other hand, is acquired from the inverted centerline velocity decay plot (Ue/Uc). Again, a line is 
extended backward from the linear self-similar decay region to the X-axis to indicate the virtual 


origin as shown in Figure 2. 
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Figure 2. Inverted Centerline Velocity Decay Plot 


B. POTENTIAL CORE REGION 


The potential core region (Figure 1), closest to the orifice, is highly influenced by the 
outlet conditions of the orifice as well as the surrounding medium. The centerline streamwise 
velocity remains approximately constant in this region; it may however exhibit a slight vena 
contracta at higher flow rates. The significant shear distortions in the flow field rapidly generate 
instabilities and create turbulent energy in the shear layers and greatly effect the initial spreading 
rate of the jet. Thought by many early investigators to be totally incoherent fluid motions, this 


turbulent flow has been proven to have ‘organized-structures’ (coherent structure) superimposed 


on its chaotic (or random) motion. This large-scale turbulence generated at the jet boundaries is 
most responsible for noise and initial mass entrainment. List/Rodi (1982) review several modes 


of wave/vortex growth and instabilities and their descriptions will not be repeated here. 


Generally, by lowering inlet core turbulent intensity (Tu) and flattening velocity profiles 
(top hat), longer potential core lengths are realized in this region. Also increasing the exit 
velocity (effectively lowering Tu and flattening the velocity profile) can enlarge the potential 


core length. 


C. INTERMEDIATE REGION 


As entrainment develops on its periphery, the potential core disappears into the 
intermediate (transitional) region. The width of this transitional region spreads rapidly outward 
as the jet progresses down stream and the average centerline velocity decreases. This re-shaping 


of the mean flow tends to force the jet into self-similarity (self-preservation). 


D. SELF-SIMILAR REGION 


Eventually, the flow attains self-similarity in which the various terms of the momentum 
equation (convection, diffusion, production, etc.) maintain constant ratios so that the various 
processes are in dynamic equilibrium, each one changing downstream at the same rate as the 
others. In this region, dimensionless velocity and length scales (a function of one geometrical 


variable only) can completely describe the flow. The flow is also said to be independent of the 


initial jet conditions including orifice shape and turbulence intensity. The width of this region 
seems to grow at a linear rate as the flow diminishes. Mean velocity profile similarity is usually 
reached within about 20-30 jet diameters downstream while turbulence velocity fluctuations 
don’t obtain “true similarity” until 70 diameters or greater depending on the type of flow. Many 
books, such as Abramovich (1963), Launder (1975) and others, discuss similarity arguments in 


much greater detail. 


Jets can obtain “exact” similarity only if they issue into stagnant or quiescent 
surroundings (as opposed to co-flow or counter-flow) and, in general, are theorized to have the 
following approximate streamwise characteristics under the assumptions of: steady flow, 


negligible body force and negligible pressure gradient (uniform static pressure). 


For axisymmetric round jets, the magnitude of axial velocity at the center of the jet (Uc) 
varies as the inverse of the axial-distance (X) from the jet exit nozzle. On the other hand, the 


radius at which the axial velocity of the jet profile is half of the centerline value (Ry) varies 


directly with X. 
Centerline Velocity Decay: Uc« >. (1) 
Half-Velocity Radius: Rin« X (2) 


Planar jets exhibit slightly different characteristics due to their two-dimensional nature. 
The centerline velocity (Uc) and transverse half-velocity length (Y1) vary in the axial direction 


(X) as follows. 


Centerline Velocity Decay: Uc« Dale (3) 


Half-Velocity Length: Yin« X (4) 


E. TEMPERATURE PROFILE 


The turbulent mixing within a jet causes a transfer of fluid particles in a direction normal 
to the main flow. The process is complex and involves momentum flowing away from the main 
stream while the surrounding fluid is being entrained toward it. In addition, there can be a 
transfer of heat or diffusion of the temperature field into the surroundings. As shown in Figure 
3, the temperature profile is considerably broader than the velocity profile due to the higher 
transverse transport of thermal energy by turbulence than momentum for gas jets. The scalar 
values of temperature, internal energy and species concentration all behave in the same way for 


these types of flow. 


Abramovich (1963), Hinze (1975), Schlichting (1979) and others describe heat transfer 
theories within the jet that are based on several simplifying assumptions. Although in general, 
the exact jet characteristics are not presented in the texts and the authors apply their (or other’s) 
experimental data to obtain the required empirical constants. However, each jet tends to be 


unique and the limited data used by each author are only applicable for that specific jet which 
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may only be similar to one or two other experiments. Even though these limited data (usually 
dated) will support the author’s general conclusions, many other experiments with comparable 
setups do not obtain the same results. Obviously not all factors (like coherent structures, etc) 
have been properly controlled in the experiments and incorporated into equations (1 to 4) and 
further investigation is certainly warranted. As will be discussed in Chapter 1, factors like 
density ratio and facility-related disturbances can greatly effect the developing characteristics of 


free jets. 


n=R(X-Xo) 


Figure 3. Typical Velocity and Temperature Profiles 
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iil. CRITICAL REVIEW OF EXPERIMENTAL RESULTS 


A. OVERVIEW 


It is essential to have reliable experimental data against which CFD predictions can be 
compared. Free turbulent jets have been measured fairly extensively throughout the years and 
plenty of data are available. The purpose of this chapter was to review these data, check its 
reliability, and select target values for a comparison with the predictions from the CFD code. As 


the reader will soon discover, this task is not as simple as it first may appear. 


B. MEASUREMENT TECHNIQUES 


In order to make judgements on the particular quality of the measurements; we will 
briefly examine the merits and shortcomings of various measurement techniques used to obtain 
the referenced data. In most free jet flows, the static pressure varies very little across the flow. 
Therefore, measurement of total pressure by a Pitot tube (or similar) can be sufficient to 
determine the mean velocity (U). The Pitot tube, due to its slow response however, can’t 
measure the fluctuating velocities that represent the most serious source of error. When 
turbulence fluctuations are high (>20% - like in the shear layer), the instantaneous yaw angles 


can be large enough to affect the probe reading without any reliable correction available. 


13 


The hot-wire anemometer has a sufficiently fast response to follow the turbulence 
fluctuations and has become the one of the more popular instruments for measuring turbulent 
flows. However, many sources of error are associated with the hot-wire measurement technique 
and not all early researchers have properly accounted for them as described by Hussein et al 
(1994) and Panchapakesan / Lumley (1993). The more common errors include improper unit 
calibration, turbulence interference between adjacent wires (multiple wire configurations), and 
failure to realize that a significant fraction of reverse flow exists in the outer shear layers of jets. 
These problems have introduced many uncertainties and the general reliability of the data 


essentially decreases with increasing turbulence level. 


The relatively new laser technologies offer much promise for high accuracy data. By 
applying a known frequency shift to the laser beams, an offset frequency is created that is 
increased or decreased by the Doppler-shift phenomena. Laser-Doppler Anemometers (LDA) 
are one of the very few velocimeters capable of measuring flow reversals. However, concerns 
still arise about proper seeding material size and dispersion required to truly represent the flow. 
Presently, only very few jet flow configurations have been re-examined using this state-of the-art 


technique and limited data are available. 


Cc. JET CHARACTERISTIC EQUATIONS 


To efficiently compare different results, the widely accepted equations given below are 
used to describe round jets in the self-similar region. Most modern researchers prefer these 


popular equation formats. Of course, some conversions were necessary to adapt all of the data 
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into this form. The velocity profiles were typically approximated using a Gaussian format 
generated from the jet’s spreading rate as shown below. The Gaussian profile tends to be a good 
approximation for the near axis regions of jets, however, the outer edges of real jets fall off 
slightly quicker as described by Malmstrom et al (1997). The following equations do not require 


that the virtual origins to be co-located. 


1. Axisymmetric Round Jets 
R, rs 
Half Velocity Spreading: fae =" 5 C (5) 
D lu D lu 
C,, = *ate (6) 
“—— D 
where: Ryu -half-velocity radius 
D -jet exit diameter 
x -axial distance from the jet exit nozzle 


Ki. —_-velocity spreading rate 
Ci, | -dimensionless velocity geometric virtual origin 
Xoju -Vvelocity geometric virtual origin 


Gaussian Velocity Profile: UL kur] (7) 
Uc 
(not all data) 
where: U -axial velocity 
Ue -centerline axial velocity 
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_ In(2) (8) 
3u 2 
K lu 
| ee eae (9) 
X-C,,D 
Centerline Velocity Decay: Ue aR x -C (10) 
Uc 2u D 2u 
es X 24 (11) 
ue D < 
where: U: -jet exit axial velocity 
Ko ~~ -centerline velocity decay rate 
Cx, -dimensionless velocity kinematic virtual origin 


Xo,2u ~-velocity kinematic virtual origin 


R 
Additional relationships: a = = (12) 
Mou lu 
R 7), 
pre ere (13) 
R i K,, 


Not all of the above variables were reported in every experiment. Sometimes the authors 
assumed the virtual origins were co-located and did not clearly indicate this within the article. 
This, in turn, constrains the reader to make certain assumptions about the jet’s characteristics, 


which may be reflected in the data presented in Tables 1 and 2. 
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As mentioned before, the temperature differential between the surroundings and the jet, 
result in temperature profiles that are wider than the velocity profiles. This same reasoning can 
also be applied to mass concentration results. Therefore, temperature and concentration are 


interchangeable and characterized in the following equations 


Half Temperature Difference Spreading: Yat _ K xX C (14) 
i It 
D 
Xx o,lt 
Ci, = 5 (15) 
where: Ry,  -half-differential temperature radius 


K;,  -differential temperature spreading rate 
Ci,  -dimensionless temperature geometric virtual origin 
Xoit | -temperature geometric virtual origin 


AT ee = ebksn’] 


Temperature Difference Profile: ATo = Te-T. (16) 
(not all data) 
In(2) 
K,. =—— 17 
3t Ke ( ) 
ee (18) 
X-C,,D 
where: T ~temperature 
T, -centerline temperature 
To -surrounding environment temperature 
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Centerline Temperature Diff. Decay: Ale. Teer. - ae (19) 
ATc Te-T,  “\D ™ 
x 20 
(om _ 0,2% ( ) 
D 
where: U. ~jet exit temperature 
Kx -centerline differential temperature decay rate 
Cz, -dimensionless temperature kinematic virtual origin 


Xo,  -temperature kinematic virtual origin 


The half-velocity (temperature) length, Rix (Rix), represents the radial distance from 
the jet’s axial centerline to the position where the velocity (temperature) is one-half the 
maximum centerline value. The values of sheets cienrae virtual origins (Cyy, Con, Cyr, C2) 
and multiplication constants (Kiy, Kou, Kit, Kz), in addition to several other relevant facts from 
several selected experiments, are listed in Table 1. As shown, the multiplication constants (K’s) 
are fairly consistent (<10% variation) between experiments. However as noted by several 
researchers, the virtual origins are found to vary widely (-0.6 < C’s < 4). Unfortunately, no 
consistent underlying reason has yet been correlated to explain the wide variations in the virtual 
origin values. In addition, not all virtual origins have been quoted in the research which makes 


comparisons between experiments difficult. 
2. Two-dimensional] Planar Jets 


The similarity analysis for the planar jet is, as expected, very similar to the axisymmetric 
round jet. However, the characteristic equations take into account Cartesian (instead of 


cylindrical) coordinates and the use of appropriate substitutions. 
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The overall resultant equation forms are: 


Y 
Half Velocity Spreading: he =k, x_ C, (21) 
H a : 
C,, == | (2) 
where: Yy, ~ -half-velocity transverse length 
H -jet height 
Gaussian Velocity Profile: UU = ob Koun] (23) 
Uc 
(not all data) 
In(2) 24 
K,, =r Cs) 
lu 
ae Y 
X-C,H | aa 
Ue Y x 
Centerline Velocity Decay: & = K,, G —~Cr, (26) 
Xx 
Cy, = (27) 
Y 
Additional relationships: yy = "4 (28) 
Mu lu 
Ym (29) 
Vy, K,, 
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Half Temperature Diff. Spreading: 


where: Yt 


Temperature Difference Profile: 


(not all data) 


Centerline Temp. Diff. Decay [ 


The values of the dimensionless virtual origins (Cyy, Cou, Ci, Co) and multiplication 
constants (Kiu, Kau, Kit, Ka) for planar jets are listed in Table 2. Again, the multiplication 


constants (K’s) are fairly consistent (<25% variation) between experiments, while the virtual 


(30) 


C = X oat (31) 


-half-differential temperature transverse length 


ia ieee lied Bee 


= =e (32) 
Ale Te-T. 
ge (33) 
Ki 
n,=|—— (34) 
a en 
Afey {Teor Y Xx 
é e=— 
eee |e =e ey eee 
or] | (5 :| es 
C,, = net (36) 
H 


origins vary widely (-6H < C’s < 4H). 
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D. VARIATION OF EXPERIMENTAL RESULTS 


Most experiments have validated the similarity analysis, as described in Chapter U, 
Section D. However, works by Kotsovinas (1976) and Jenkins / Goldschmidt (1973) have 
continued to shed doubts on the linear behavior of planar jets in the far similarity region. 
Depending on the downstream range investigated, the planar jet centerline decay data appears to 
be slightly non-linear (greater) when using the standard similarity equations (26 & 35) far from 
the jet exit (X/H > 200). Therefore, the regression range used by researchers to obtain virtual 


origins and multiplication constants are important. 


1. Incomplete Data 


Accurate universal jet predictions are very hard to acquire. Several factors can influence 
the development of each jet, as well as describe their general shape. Almost no text or 
experimental report has incorporated all of the required factors to completely characterize jet 
flow. Each is based on very limited data and a general overall behavior is assumed for all 
comparable jets. Although this is partially true, due to the similarity relationships, one can not 


apply these results to all jets (even when flow conditions seem to be the same). 


The most commonly recognized difference among experiments is the location of the 
virtual origins. Many researchers have not clearly recognized that in some situations the 
geometric (Cjy, Cy) and kinematic (Coy, C2) virtual origins do not necessarily coincide as would 


be expected by the similarity arguments [See Goldschmidt / Bradshaw (1981) and Flora / 
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Goldschmidt (1969)]. This is further complicated when some authors fail to report all of the 
virtual origins. This consensus usually arises from the assessment that only the multiplication 
constants (K’s) are required to adequately describe jet characteristics. However, one can not 
completely describe the longitudinal position of the centerline decay and lateral spread profiles 


without knowing the reference points (virtual origins) upon which they’ re based. 


The virtual origin, therefore, is an important parameter, although it may be quite 
challenging to predict. Papadopoulos / Pitts (1999) have developed successful centerline 
prediction equations for iso-thermal jets using the initial momentum flux and local centerline 


turbulence intensity. However, further research on heated jets has not been performed. 


2. Jet Variations 


To assume that the results of one experiment will completely match another is erroneous. 
Slight variations in the measured multiplication constants and virtual origin values (See Table 1 
and 2) can have significant effects on jet predictions using the similarity equations (Equations 5- 


36). 


In addition, free jets have been found to contain large-scale coherent disturbances within 
the shear layers on top of the seemingly random turbulent motions [See Thomas / Prakash 
(1991), Moum ef al (1983), Hussain (1983), Gutmark / Ho (1983) and Crow / Champagne 
(1971)] which are initial condition dependent. These disturbances can be excited relatively 


easily and can drastically effect the entire motion of the jet (non-quasi-steady). Hussein et al 


22 


(1994) and George (1989) further described how a more general analysis of the similarity theory 
reveals that a self-preserving state can, indeed, be uniquely determined by the initial conditions. 
Hence, a family of farfield similarity solutions could exist for each type of jet, depending on its 


initial conditions. 


The density (and/or temperature) differential effect on jet development has been 
investigated by Drobniak et al (1998), Sautet / Stepowski (1995), and Russ / Strykowski (1993). 
Further articles by Kyle / Sreenivasan (1993) and Monkewitz et al (1990) describe how “hot jet 
to surrounding environment” density ratios less than approximately pje/p. < 0.6 can cause a self- 
excited large-scale coherent mode that produces a marked jump in the spreading rate. This mode 
is caused by the formation of a very regular sequence of intense vortex ring structures within the 
jet. The behavior of the mode has been found to be independent of background disturbances 
within typical experimental Reynolds Number ranges (10”2 to 10S). The resulting “strong 
pairing process” of the mode leads to abnormally large centerline velocity fluctuations and early 
potential core breakup. Current commercial CFD codes, solving RANS equations, cannot model 
these excited modes which have become major obstacles when accurately predicting free jet 


flows. 


Some additional experimental variables that haven't always been controlled or reported 
include the following. Fan vibrations felt at the nozzle and surrounding area disturbances are 
known to cause significant changes in the jets, especially at low flow rates as described by 


Gutmark / Ho (1983). Modern researchers have tried to prevent such interactions to varying 
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degrees of success, but their reassessment of earlier experiments are still quite limited. Small or 
semi-confined test sections and inadequate equipment have also tainted some earlier experiments 
as described by Malmstrom et al (1997), Hussein et al (1994), and Panchapakesan (1993). The 
overall room size and re-circulation patterns seem to deform and shift the immerging jet enough 
to effect the results between different experimental facilities. Extreme care has been taken by 


modern researchers to minimize these effects. 


24 


Sé 


Table 1. Axisymmetrical Round Jet Experimental Data 


Investigator Jet Test Pd Po Jet Jet 
Dia Range Ratlo Vel. Tu 
(m) Prof. 
(D} {XD} {%} 


Hussein, Capp & 0.0254 | LDA | 10-120 

George (1994) 

Panchapakesan & 0.0061 | FHW Top 

Lunley (1993) Hat 
Wygnanski & 0.0264 Hot 0 Top 
Fielder (1969) Wire Hat 


Corrsin & Kistler 0.0127 20-50 73x10 ue some 0.0945 
1954 Hat 


TEMPERATURE OR CONCENTRATION DIFFERENCES 
Grandmais, 0.071 Laser | 10-80 27410 Cone. . : 
Rathgeb. & Becker Tech. Effect 
1982 
List & Rodi Top 0.089 87.3 0.105 0.194 
(1982) Hat 
Data Review 


Notes: t: (14-12.127' £28159 bi] H: (14.63.78?) ” tH: ( U )~ T/C: thermocouple FHW: Flying hot wire 
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Table 2. Two-dimensional Planar Jet Experimental Data 


Investigator Test PelPo Jet Jet . Spread 
Range Ratlo Vel. Tu : i i Virt. 
i Prof. . é Orig. 
{H} {X/H} {Cry 


ISOTHERMAL CONDITIONS 


Thomas & Pradkash | 0.0127 Hot 0.8x10 Top 0.24 0.14 0.22 Pes 
(1991) Wire Hat 


Goldschmidt & 0.0098 Hot 10-80 2.7x10" F 
Bradshaw Wire . , ‘ , é 
(1981) Top , . F I i 

Hat : : 
Hussain & Clark 0.0318 Hot Top | 03 0.113 -2.11 
(1977) Wire Hat 
Gutmark & 0.0127 Hot i neil Top | 0.2 0.188 
Wygnanski(1976) Wire Hat 
Flora & 0.0127 2.6 x10" Top | 1.06 0.130 -15.0 0.227 2.00 
Goldschmidt (1969) Hat 1.2) 0.113 6 0.144 y 

1.28 0.127 -5.00 0.222 

Miller & Comings 0.0127 Hot Top 0.097 -1.57 167 5 AT 
(1957) Wire Hat 
van de Hegge 0.005 1 eae Top 0.096 0.163 -0.6 75 
Zijnen (1957) 0.010 Hat 0.099 0.205 -17 70.7 


TEMPERATURE OR CONCENTRATION DIFFERENCES 


Chen & Rodi 
(1980) Data Review 


Davies, Keffer & Ia 
Baines (1975) Cold 
Wire 
Jenkins & 0.013 Wire | 10-55 1.4x10' 0.96 0.088 
Goldschmidt (1973) and 0.94 0.091 
Pitot 0.90 0.096 
van der Hegge 0.005 Pitot | 2-40 1.3x10 0.84 
ear 
TiC 


Notes, +: AGARD Table t: “ we -x,(4-c,) ATe ttt: AZ _ § 430m? +22007,4 -300007, pF” 


E. OVERALL TRENDS 


- Several overall trends in experimental data have been noticed throughout the years with 
varying degrees of consistency. Although a particular reported trend might have an exception in 
one experiment, the overwhelming inclinations (as described by many additional researchers) 


will be described here. 


1. Exit Velocity Effects 


Several earlier experiments report a diameter-based Reynolds Number (Reg #) 
dependence of jet flows with exit Reg #’s less than approximately 5x10*. By varying the jet exit 
velocity, Reg # dependence was seen as a change in the spreading rate and centerline decay. A 
more comprehensive study by Malmstrom et al (1997) which involved varying the exit velocity, 
as well as the jet diameter, concluded that the variance was based on jet exit velocity, not the Reg 
#. As jet exit velocity was increased beyond 6 m/s, the velocity centerline decay multiplication 
factor (K2,) decreased to an asymptotic 0.17 value (for round jets) which agrees well with Table 
1. Below 6 m/s, as jet exit velocity is decreased, the centerline decay and spreading rate 


increase. 


2. Turbulence Intensity Effects 


The dominant factor that seems to effect experimental virtual origin location is turbulence 
intensity at the jet exit profile. Both the centerline turbulence intensity and boundary layer 


thickness can drastically shift the origins in seemingly random directions. Several researchers 


2) 


have attempted to qualify these effects with limited success in narrowly defined jet 
configurations. Changes in the turbulence intensity have also sometimes masked the effects of 
velocity dependence on jet development. As the jet exit velocity increases, the core turbulence 
intensity generally drops (unless artificially held constant) as the boundary layers become 
thinner. Therefore, as exit turbulence intensity decreases, the jet spreads less and the centerline 
decay decreases. However, Goldschmidt / Bradshaw (1981), Kleis / Floss (1974), and Flora / 
Goldschmidt (1969) have shown that the above Tu effects still hold true when the exit velocity is 


held constant. 


F. TARGET VALUES 


In order to compare the CFD results to experiments, target values must be chosen. A 
band (or range) of experimental data was select over a specific experiment because every jet 
tends to be unique as described in Section (D). After analyzing Tables 1 and 2, evaluating the 
references, and studying more recent investigations, the following parameters were chosen to 


compare and model typical jet flow for this research. 


1. Centerline Decay and Spreading Values 


As shown in Tables 1 and 2, the spreading rates and centerline decay values are fairly 
constant, but do exhibit a slight variation (~10% for the round jet and ~25% for the planar jet). 
Based on measurement technique, data scatter, test range, researcher credibility, the following 


values will be used to compare CFD results. 
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Table 3. CFD Simulation Target Values 

Temp | Spread 
Virt. 
Orig. 


Vel 
Spread 
Rate 


Virt. 
Orig. 


Flow Type 


Orig. 
Kau Cy Koy Cu Kir Cr Ka Cr 
Round Jet | 0.085- | -1.0to | 0.160- | -1.0to | 0.098- 0.0 to 0.185- 0.0 to 
0.095 4.0 0.172 4.0 0.110 2.0 0.194 2.0 
Planar Jet | 0.096- | -6.0to | 0.155- | -2.0to | 0.125- | -5.0to 0.250- | -5.0 to 
0.120 2 0.220 40 0.140 0.0 0.280 0.0 
6H 5H 


Pies [sa [ie |S _ 


2.0 
[Var | 25% | 8H [| 419% | CH 


Target virtual origin ranges were selected to include the extreme reported variations. As 
mentioned, current steady-state Favre-Averaged Navier-Stokes (FANS) CFD codes can not 
model many important external factors (i.e., nozzle vibration, room disturbances, etc.) that can 
significantly effect the virtual origin position. Therefore, the siiial origin target range is used 


only for vague simulation result comparisons. 
2. Turbulent Prandtl Number 


The turbulent Prandt! Number (Pr,) quantifies the transport rate of momentum to that of 
heat and was used as an input to the CFD code k-e turbulence / energy model. Free jet flows, 
which are not effected by walls, display different Pr, variations normal to the primary flow 
direction than typical internal flows (i.e. equilibrium pipe flow). The turbulent Prandtl Number 
varies across the round jet’s profile and was quantified by Reynolds (1976) as shown in Table 4 


and Figure 4. 
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Table 4. Turbulent Prandtl Numbers for Jets 


Local Pr; Recommended 


Global Pr; 


Flow Type 
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Figure 4. Typical Prt and y Variations across a round and planar jets 
[Based on Reynolds (1976)] 
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The Turbulence intermittency factor (y) in the jet profile varies from unity near the jet 
centerline (¢) and inner part of its shear layer to 0.5 and below in the outer edges of the shear 


layer. The most striking feature of the results is the difference between the Pr, in the body (y=1) 
and the outer edges of the highly turbulent shear layer. A Pr greater than unity indicates that 
momentum is being transferred more readily than heat in these highly intermittent regions via 


relatively large-scale turbulent structures. 


These variations, although significant across the whole jet, are fairly small across the 
largest inner portion of the jet and can be approximated using a single global value for the entire 
jet. In general, the body values for round jets are larger than for planar jets. Chen / Rodi (1980) 
and Reynolds (1976) recommend the well accepted global Pr,‘s listed in Table 4 when the 
density differences are small and the molecular Prandtl Number (Pr) is close to unity (Priair ~ 
0.7). These global approximations, of coarse, result in small CFD temperature over-prediction 


> 


errors in the outer edges of the jets. 
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IV. CFD SIMULATION 


A. CODE AND COMPUTER 


The CFD Research Corporation (CFDRC) developed the computational fluid dynamics 
package used in this research to analyze the different jet configurations. The CFD-ACE+ 
program (version 6.2) was run on an IBM/PC-based Gateway (Model GP6-450) 450 MHz clock- 
speed computer. The Random Access Memory was upgraded to 384MB to handle grid files with 
up to 300,000 three-dimensional unstructured cells. A brief description of the CFDRC code and 


its operation is provided in Appendix A. 


B. THE PROCESS 


To gain experience with the code, isothermal jets were initially modeled using small 
domains with a coarse grid to minimize solution run times. This process allowed the flow 
conditions to be fully understood before the complexity of heat transfer was added. After the 
code was mastered at each step, an additional feature.or component was activated. Grid 
refinement was also developed in the high gradient areas. Only the final heated jet models are | 
discussed within this report. In the anticipation of further research with more complex flows 
using this code, only 3-dimensional unstructured tetrahedral cells were used in the generation of 


all grid domains. 
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C. HEATED AXISYMMETRICAL JET MODEL 


The heated axisymmetrical round jet was modeled to investigate the ability of the code to 
predict these types of flows. Although this code can not model all factors effecting jet 
development (i.e., coherent structural modes, etc.), an overall performance evaluation of the code 


compared to experimental data was desired. 


1. Grid Domain 

A one-quarter section axisymmetrical grid was used to accurately model the developing 
jet and its surroundings while minimizing the total grid size. Smaller sections (i.e. 1/6, 1/8, etc.) 
were not used to prevent possible symmetric edge interference problems. Typical grid 
dimensions and overall configuration are shown in Figure 5. The total grid length was selected 
to be 40 jet exit diameters with a grid radius of slightly less than 8 jet exit diameters in order to 
capture the flow features without boundary interference. The overall grid domain was divided 
into two separate sub-domains. The inner sub-domain consisted of smaller cells to more 
accurately capture the high gradients of the main jet flow and its free shear layer (adaptive 
meshing was not available). The inner sub-domain was constructed as an inverted truncated 
cone to mimic the jet’s radial spread as it travels downstream. The inner sub-domain 
encompasses most of the jet’s flow throughout the entire grid domain. The jet enters this sub- 
domain through the lower face and exits through the top with entrainment coming in from the 
sides. The bottom boundary of the inner sub-domain is slightly larger than the jet exit diameter 


to facilitate a finer mesh (more accurate solution) around this transitional flow area. 
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__ Jet Exit 
D/2=0.0127m 


, 


Figure 5. Axisymmetrical Jet Grid 


Before generating the internal tetrahedral volume grid, the inner sub-domain’s boundary 
surface grids were automatically generated using the following dimensionless control 


parameters. 


* Maximum cell size: 0.007 
* Minimum cell size: 1.143x10° 
« Curvature Criterion: 8 


« Surface Transition Factor: 1.1 
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The Curvature Criterion defines the approximate maximum angle of curvature per cell 
that will be generated at edges that are curved. The Transition Factor limits the rate of cell 
growth (size change) when generating triangular surface cells from a curved edge to areas far 


from all curved edges. 


After all surfaces of the inner sub-domain were gridded, a separate volume tool was used 
to invoke the tetrahedral volume grid generation algorithm. The tetrahedral volume grid 
generation algorithm used in CFD-GEOM is a variation of the advancing front grid generation 
technique. The size of the generated tetrahedrons change smoothly, based on the properties of 
the surface grids. One additional mechanism (besides those used in surface generation) was used 
to contro] the tetrahedral grid volume generation. The Volume Transition Factor, similar to the 
Surface Transition Factor, determined how far off a face that a field point would be generated for 
the formation of a new tetrahedron. For the inner sub-domain, the Volume Transition Factor was 


set at 1.1 


In general, the outer sub-domain contained larger cells to model the airflow in the 
relatively slower moving surrounding environment. The outer boundaries of this domain were 
made vertical on the sides and horizontal on the top and bottom to simplify boundary conditions 


at these areas. The grid generation controls for this domain were set as follows. 
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« Maximum cell size: 0.05 

* Minimum cell size: 1.143x10° 
¢ Curvature Criterion: 30 

¢ Surface Transition Factor: 1.2 


¢ Volume Transition Factor: 1.2 


The outer sub-domain grid was generated after the inner sub-domain so that the 
tetrahedral cells could progress smoothly from the small inner cells to the larger outer sub- 


domain cells at the grid boundaries. 


Overall, the entire domain contained 317,273 tetrahedral cells with the following grid 
qualities as computed by the code. The CFDRC User Manual (2000) can be consulted for 


further definition of each criterion. 


+ Dihedral Angle: 3 cells were less than 5 degrees (default min.=5 degrees) 
° Worst Centroid/Face Angle: 7.24 degrees (default min.=5 degrees) 
¢ Sliver Quality: 5 cells were greater than 7 (default max.=7) 


* Skew Quality: 52 cells were greater than 0.95 (default max.=0.95) 
Only a small faction of the cells in non-critical areas violated the default limits and the 


gridded domain was determined to be satisfactory. Several other complex domain configurations 


were tried to further concentrate small cells in high gradient areas without greater grid quality 


af 


success. In addition, the higher domain complexity adversely effected the solution convergence. 
More controls are needed to properly manage the automatic grid generation process for 


unstructured cells and CFDRC plans to release future versions addressing this problem. 


2. Solver Setup 

a. Overall 

Table 5 below provides a summary of all solver [CFD-GUI / CFD-ACE (U)] 
inputs. To simulate the heated jet, the flow, heat transfer and turbulence modules were activated. 
The relatively simple k-e turbulence model was chosen to gain experience manipulating the code 
while maximizing solution convergence. As mentioned in Appendix B, the k-€ turbulence model 
constants can not be changed (within the code) for axisymmetrical flows as suggested by 
Launder/Spalding (1974). Therefore, without being able to vary Cy and C,; (turbulence model 
closure constants) across the profile, the overall prediction accuracy for axisymmetrical jets was 


expected to be poor. 


b. Volume Conditions 

Air at ambient conditions was selected as the computational fluid. To simulate 
the heated jet (AT = 30K), density calculations evaluated using the Ideal Gas Law while the 
viscosity was determined by Sutherland’s Law [U=AT*7/(B+T)], where A and B are constants 
specified in Table 5. Since the simulation involved small differential temperatures, fluid 
conductivity (K) and specific heat (Cp) evaluations were determined by selecting the constant 


molecular Prandtl Number (Pr = 0.7) option and constant Cp = 1000 J/kg-K, respectively. 
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C. Boundary Conditions 

The grid inlet boundary was modeled as 0.0254m (1-inch) diameter jet exit nozzle 
with a constant velocity (top hat) profile. Since the grid domain was modeled as a 1/4" section, 
the inlet boundary is shaped like a 1/4" “pie section” with a radius D/2 = 0.0127m. The ideal 
‘top hat’ outlet jet velocity profile approximates the flow that most researches try to develop 
using large nozzle contractions. Although ‘real’ jets have small boundary layers which can 
drastically effect their development, an ideal top hat profile was used to simplify the CFD model. 
With a properly supplied sub-routine, a user can define any desirable inlet jet velocity profiles, 
however, the improvement in the final CFD solution was expected to be negligible when 


compared to the widely varying experimental jet data. 


The constant jet nozzle exit velocity (U) {grid inlet} was initially set at 40m/s, in 
the axial direction (x) , with a uniform temperature of 330 degrees, K. All inlet turbulence 
kinetic energies (k) needed for the Favre-Averaged Navier-Stokes (FANS) conservation 
equations (Equations 42 and 43) were calculated from Equation 37, in which the turbulence 
intensity (T,) was assumed to be contemporary experimental values. A turbulence intensity of T, 


= 0.5% was originally selected which resulted in k = 0.06 m/s’. 


k= =U? (37) 
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The inlet dissipation rate (€) was input into the code by using the ‘length scale 
option’ as discussed in Appendix B. This method required that the inlet size dimension (L) be 
specified, then the turbulent length scale (0.3% of the size of the inlet size dimension) is 
subsequently calculated within the program. Since the grid inlet boundary is a “pie section”, the 


size dimension was set at L = 0.0127m (inlet boundary radius). 


The straight sides of the 1/4 section grid domain were set as symmetric 
boundaries. The rest of the grid boundaries (top, curved front surface and rest of the bottom) 
were established as fixed (static) ambient pressure outlets. This setup, however, does allow flow 
to enter the grid domain through the side and bottom boundaries for entrainment. The outlet 
boundary parameters were set to characterize the flow as it entered the grid domain from the 


infinite surroundings. 


d. Initial Conditions 
The initial conditions are displayed in Table 5. The starting turbulent viscosity 
was artificially elevated (via k and L) to enhance solution convergence during the earliest 


solution iterations. 


é. Solver Control 
A 1*-order upwind spatial differencing scheme was chosen for all variables in the 
jet models to simplify convergence operations and minimize solution times (typically 2-3 days). 


The default CGS+Pre (conjugate-gradient-squared plus preconditioning) linear equation solver 
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was also chosen to resolve the algebraic equations for each dependent variable and the 
controlling parameters. The solver performed satisfactory without approaching its internal 
default maximum number of sweeps criterion (50 or 500 sweeps, depending on the variable). 


The under-relaxation parameters were also left at their default values (see Table 5), which 


seemed to work adequately. 
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Table 5. Summary of Round Jet Solver Settings 
Problem Type 


Flow 
Heat Transfer 
Turbulence 


Model Options 


Domains Density Evaluation Method Ideal Gas Law 
M.W. = 29 (default 


Viscosity Evaluation Method Sutherland’s Law 
A= 1.4605E-6 
B= 112 (default 
Specific Heat Evaluation Method Constant 
Cp = 1000 J/kg-K (default 
Pr=0.7 
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Under-Relaxation Parameters 


lEnthalpy 0.05 (default 
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0.2 (default 
[Pressure | efauit) 
1 (default 
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D. HEATED PLANAR JET MODEL 


1. Grid Domain 

Again, only 3-dimensional unstructured tetrahedral cells were used in the generation of 
the planar grid domain. A rectangular grid was used to model half (1/2) of the developing planar 
jet and its still surroundings. See Figure 6 below for typical grid dimensions and overall 
configuration. The total grid length was selected to be 40 jet exit heights (H). A grid width of 
approximate 11.5 jet exit heights and depth of 4 jet exit heights was established in order to 
capture the flow features without boundary interference. The overall grid domain was again 
divided into two separate sub-domains. The inner sub-domain consisted of smaller cells to more 
accurately capture the high gradients of the main jet flow and its free shear layer. The jet enters 
the inner sub-domain through the lower face and exits through the top with entrainment coming 
from the non-symmetrical far side. The bottom boundary of the inner sub-domain is slightly 
larger than the jet exit diameter to facilitate a finer mesh (more accurate solution) around this 


transitional flow area. 
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Jet Exit 
(H/2) = 0.0065m 


Figure 6. Two-Dimensional Jet Grid 


The inner sub-domain’s boundary surface and volume grid was automatically generated 


using the following control parameters. 


* Maximum cell size: 0.004 
* Minimum cell size: 5.2x10° 
* Surface Transition Factor: 1.1 


* Volume Transition Factor: 1.1 


The outer sub-domain grid generation controls were as follows. 


¢ Maximum cell size: 0.01 
¢ Minimum cell size: 5.2x10° 
¢ Surface Transition Factor: 1.1 


¢ Volume Transition Factor: 1.1 


The outer sub-domain grid was generated after the inner sub-domain so that the 
tetrahedral cells could progress smoothly from the small inner cells to the larger outer sub- 
domain cells at the grid boundaries. The typical cell size was smaller than the axisymmetric jet 


model due to the smaller overall grid volume (limited to approximately 300,000 cells). 


Overall, the entire domain contained 297,657 tetrahedral cells with the below grid 
qualities as computed by the code. All grid quality checks were within limits and the gridded 


domain was determined to be satisfactory. 


¢ Worst Dihedral Angle: 6.54 degrees (default min.=5 degrees) 
* Worst Centroid/Face Angle: 18.11 degrees (default min.=5 degrees) 
¢ Sliver Quality: 1 cell at 7.01 (default max.=7) 


* Worse Skew Quality: 0.93 (default max.=0.95) 
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2. Solver Setup 
a. Overall 
Table 6 below provides a summary all of the solver inputs, which are very similar 


to the round jet simulation. 


b. Boundary Conditions 

The 2-dimensional jet exit nozzle (grid inlet boundary) was modeled with a height 
of 0.013m (1/2-inch) and a constant ‘top hat’ velocity profile. Since the grid domain was a 1/2 
section, the inlet boundary has a half-height of D/2 = 0.0065m. Therefore, the inlet dissipation 
rate (€), using the length scale option, was set at L = 0.0065m to match the grid inlet size (similar 


to the axisymmetrical jet model). 


The bottom boundary was established as an adiabatic wall to emulate the 
configuration of most heated planar jet experiments. The rest of the grid boundaries were 
established as fixed ambient pressure outlets. The outlet boundary parameters were again set to 


characterize the surrounding flow if it entered the grid domain from the infinite surroundings. 


c. Initial Conditions 
The initial conditions are displayed in Table 6. The starting turbulent viscosity 
was artificially elevated (via k and L) to enhance solution convergence during the earliest 


iterations. 
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d. Solver Control 

A 1“-order upwind spatial differencing scheme was chosen for all variables in the 
jet models to simplify convergence operations and minimize solution times (typically 3 days). 
The default CGS+Pre (conjugate-gradient-squared plus preconditioning) linear equation solver 
was also chosen to resolve the algebraic equations for each dependent variable and the 
controlling parameters. The solver performed satisfactory without approaching its internal 
default maximum number of sweeps criterion (50 or 500 sweeps, depending on the variable). 
Again, the under-relaxation parameters were also left at their default values, which seemed to 


work adequately. 
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Table 6. Summary of Planar Jet Solver Settings 


Problem Type 


Flow 
Heat Transfer 
Turbulence 
Global Steady Time Dependence —~—“C~s—SCSCsd 
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Solver Control 


Spatial Differencing Scheme Upwind (default 


Density Upwind (detautt 
Enthalp Upwind (default 
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CGS+Pre (default 
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V. COMPARISON OF CFD RESULTS TO EXPERIMENTS 


CFD simulation velocity and temperature profiles at different axial distances for both jets 
are displayed in Figures 7-8 for the round jet and Figures 15-16 for the planar jet. By using the 
similarity equations (Equations 5 through 36), the non-dimensional plots (Figures 9-14 for round 
jet and Figures 17-22 for planar jet) were generated to investigate the self-similar ranges of each 
jet. As shown, the simulated round jet doesn’t achieve self-similarity until approximately 30 
diameters downstream of the jet exit (see Figures 9, 10, 13, and 14). The simulated planar jet, on 
the other hand, reached similarity within approximately 25 nozzle heights (see Figures 17, 18, 


21, and 22). 


The centerline decay and spreading rates, as well as the virtual origins, for planar and 
round jet simulations are listed in Table 7. Regression ranges were calculated in the similarity 
region between 30 and 40 jet diameters (or slot heights) for both jet types. Several values (half- 
radius/length or centerline) were analyzed by a Microsoft Excel spreadsheet regression program 
to determine the slope [spreading or centerline decay rate (K)] and imaginary X-intercept [virtual 
origin (C)] of each plot. See Appendix C for a typical sample of the Excel spreadsheet used for 


round jet calculations. 
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Figure 7. Round Jet Velocity Profiles 
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Figure 8. Round Jet Temperature Profiles 
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Figure 9. Round Jet Centerline Velocity Decay 
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Figure 10. Round Jet Centerline Temperature Decay 
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Figure 11. Round Jet Velocity Profile in the Self-similar Region 
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Figure 12. Round Jet Temperature Profile in the Self-similar Region 
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Figure 15. Planar Jet Velocity Profiles 
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Figure 16. Planar Jet Temperature Profiles 
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Figure 17. Planar Jet Centerline Velocity Decay’ 


CFD [K2t=0.209, C2t=9.1] 
~ ——-High Target [K2t=0.280, C2t=-5] 


Figure 18. Planar Jet Centerline Temperature Decay 
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Figure 19. Planar Jet Velocity Profile in the Self-similar Region 


CFD [K1t=.149, C 1t=8.6] 
High Target [K1t=0.140, C1u=-5} 
€ 
r 
o 
Ee 
= 
Fe 
c 
0 0.05 Of O15 O2 025 O03 0.35 
Eta=Y/(X-C1tH) 
Figure 20. Planar Jet Temperature Profile in the Self-similar Region 
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Table 7. CFD Simulation Predictions 


Spread 
Virt. 


Orig. 


Round So 
Values 0.095 4.0 0.172 4.0 0,110 2.0 0.194 2.0 
Planar | CFD | 0.092 | 3.7 | 0.154 | 52 | 0.149 | 86 | 0.209 | 9.1 | 


Jet Target | 0.096- | -6.0to | 0.155- | -2.0to | 0.125- | -5.0to | 0.250- | -5.0 to 
Values | 0.120 2.0 0.220 4.0 0.140 0.0 0.280 0.0 


The CFD predictions can be compared to the target values listed in Table 7 and 
visualized in Figures 9 through 12 for the round jet and 17 though 20 for the planar jet. These 
non-dimensionalized CFD profile plots were generated at X/D=40 for the round jet (X/H=35 for 
the planar jet). The target value ranges, indicated on the figures for each jet, represent the 
reported variation of experimental data (target value range) in the self-similar region only. The 
‘High Target’ curves indicate the upper bound on experimental data while the ‘Low Target’ 
curves represent the lower experimental data boundaries. As indicated, none of the CFD 
predictions fell within the targeted ranges. Evaluations were performed to quantify the 
difference between CFD results and the target value ranges. Each CFD value (K or C) was 
compared to the ‘middle value’ of the target value range and presented as the ‘Eval’ values in 


Table 7. 
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A. HEATED AXISYMMETRICAL JET RESULTS 


As expected, the k-e turbulence model simulated the round jet poorly as discussed in 
Section IV.C.2 and Appendix B. The round jet results indicated that the code over-predicted the 
velocity-spreading rate by 24% and the temperature-spreading rate by 29% when compared to 
the middle of their respective target value range. In addition, the centerline velocity and 
temperature decay rates were also over-predicted by 21% and 30%, respectively. The geometric 
and kinematic virtual origins were over-predicted, as well, by approximately 75 diameters for 
the velocity profiles and 10.5 diameters for the temperature profiles. Although experimental 
virtual origin data display significant variation, typical values are usually reported close to the jet 


exit [-1 < X/D <4]. 


B. HEATED TWO-DIMENSIONAL JET RESULTS 


Planar jet simulation was generally closer to experimental data ranges, with an under- 
prediction of the velocity-spreading rate of approximately 17% with an over-predicted 
temperature-spreading rate of 12%. The centerline velocity and temperature decay rates were 
both under-predicted at 22% and 27%, respectively. Again, the geometric and kinematic virtual 
origins were over-predicted by approximately 7.5 slot heights for the velocity profiles and 10.5 


slot heights for the temperature profiles. 
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C. DIMENSIONAL JET COMPARISONS 


To evaluated the actual differences between the CFD results and experiments, the non- 
dimensional experimental parameters (K’s and C’s) were used to generate actual profiles 
[velocity (m/s) and temperature (K: degrees Kelvin)] using the similarity equations (Equations 5 
through 36). By comparing the actual profiles, further insight was discovered about the accuracy 


of the CFD code. 


The ‘middle values’ of the target value ranges (experimental data ranges) were chosen for 
comparison with the CFD profiles (see Table 8). These middle values were used to generate 
‘average experimental profiles’ and were not based on any specific study, but did establish a 


foundation for evaluation. 


Table 8. Simulation Comparison Parameters 


Spread 


Round 
Jet 


Ki Cw Kou Cry Kit Cr Cx 
cel cl Dl Ell aad lied al 
Target 


0.092 0.154 0.149 | 8.6 | 0.209 


Target 
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The CFD jet simulation exit velocities [Ue=40m/s round jet (Ue=35m/s planar jet)], 
temperatures (Te=330K), and jet sizes [D=0.254m round jet (H=0.013 planar jet)] were used as 
inputs to the produce the average experimental profiles. This approach, of coarse, ignores the 
effects of velocity and temperature on the development of coherent structures within free jets. 


However, an overall comparison was still desired. 


The CED and average experimental profiles were evaluated at X/D = 40 for the round jet 
and are displayed in Figures 23 and 24. The higher CFD centerline values (Uc, Tc) are the result 
of lower Ue/Uc and ATe/ATc values at X/D = 40 as displayed in Hanes 9 and 10. Even though, 
the CFD results displayed higher centerline decay rates (Kou, Kz), their excessively large virtual 


origins (C2, C2) produced higher centerline values (Uc, Te) at X/D = 40. 


The planar jet profiles were evaluated at X/H = 35 and are shown in Figures 25 and 26. 
Again, the CFD simulations display larger centerline values (Uc, Tc) due to excessive virtual 
origins (see Figures 17 and 18). Even though the planar jet simulation predicted spreading and 
centerline decay rates (K’s) closer to experimental data than the round jet simulation, the 
differences in the actual profiles at X/H = 35 was generally worse for the planar jet due to its 
excessive virtual origins. At different distances from both jet exits, the actual results will vary 


due to the offsetting errors in the virtual origins and spreading/decay rates. 
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Figure 23. Actual Round Jet Velocity Profiles at X=1.016m (X/D=40) 
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Figure 24. Actual Round Jet Temperature Profiles at X=1.016m (X/D=40) 
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Figure 25. Actual Planar Jet Velocity Profiles at X=0.455m (X/H=35) 
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Figure 26. Actual Planar Jet Temperature Profiles at X=0.455m (X/H=35) 
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VI. SENSITIVITY ANALYSIS 
To evaluate the sensitivity of the CFD models and possibly produce more accurate 
results, various parameters were changed to observe their resultant effects. This baseline 


knowledge will be used in further research to help predict jet-type flows more accurately. | 


A. JET TURBULENCE INTENSITY 


The jet turbulence intensity was varied from 0.002 to 0.01 on the round jet (nominal = 
0.005) and from 0.001 to 0.01 on the planar jet (nominal = 0.002). Table 9 and Figures 27 to 34 
in Appendix D (Section A) display the results. The round jets virtual origins (C’s) decrease 
slightly with increasing jet turbulence intensity while the multiplication constants (K’s) remain 
relatively constant. The planar jet displays similar characteristics, but deviated in different 
directions. The smaller differences may be due to slight variations in the regression evaluation 
and not a result of jet turbulence intensity changes. Overall, the effects were minimal (<5%) 


within the given turbulence range. 


B. TURBULENT PRANDTL NUMBER 


The turbulent Prandtl Number (Pr,) had one of the largest effects (~15%) on the jet 
predictions. The effects of changing the Pr, by +0.2 are shown in Table 10 and Figures 35 to 42 
in Appendix D (Section B). For both jets, the Pr, had no effect on the velocity characteristics as 
expected, but markedly changed the temperature features. When the Pr, was lowered, the CFD 


code correctly produced wider temperature profiles (increased spreading) and larger centerline 
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temperatures decay rates as more emphasis was artificially placed on the heat transport 


mechanism relative to the momentum transport mechanism. 


In general, the Pr, can be adjusted to produce the desired temperature characteristics in 
relation to the velocity profiles. For both jet simulations, the recommended free jet Pr,’s (0.5- 
Planar, 0.7-Round) produced temperature profiles that were approximately same magnitude 
greater than the velocities profiles as reported by experimental data. However, for the round jet 
especially, the predicted velocity profiles were excessively wide, therefore the resultant 


temperature profiles were also too broad. 


C. JET LENGTH SCALE 


Varying the characteristic length (L) artificially changed the jet turbulent length scale (£) 
as described in Appendix B. The characteristic length was varied from 0.001m to 0.1m for both 
jets to investigate the resultant effect. Table 11 and Figures 43 to 50 in Appendix D (Section C) 
display the results. As shown, changes in the results were negligible (<1%) and the original 


length scales used in each simulation were determined to be adequate. 


D. JET VELOCITY 


The jet velocity was varied from 20 to 60 m/s on the round jet (nominal = 40m/s) and 
between 15 and 55 m/s on the planar jet (nominal = 35m/s). These ranges were limited to typical 


experimental values in order to identify any trends at these lower velocities. The results are 
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shown in Table 12 and Figures 51 through 58 in Appendix D (Section D). The change in 
velocity did not appreciably effect (<1%) either jet prediction. Since no boundary layers 
alterations were simulated at the jet nozzle (top hat velocity profile), the velocity variations did 


not change the solution results. 


E. JET TEMPERATURE 


The jet exit temperatures for both models were elevated to 500K degrees (200K above 
ambient surroundings) to show the effects of a larger temperature difference. By limiting 
differential temperatures to less than 200K, errors associated with using the “constant specific 
heat capacity” option within the simulation were minimized. As mentioned before, the “constant 
specific heat capacity” option greatly reduced solution convergence problems and is consistent 


with low working differential temperatures. 


The results are displayed in Table 13 and Figures 59 to 66 in Appendix D (Section E). 
As indicated the temperature difference effects on the velocity and temperature profiles (and 
spreading rate) are negligible (<1%) within the simulated temperature range. However, the 
velocity and temperature centerline decay rates (Kz, Kz) increased markedly (~5%). The 
centerline decay characteristics are consistent with experimental data as described in Chapter II, 
Section C.2, however the simulated spreading rates (Ki, K;,) do not exhibit the coherent 


structure effect, as expected. 
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F. GRID REFINEMENT 


The individual cell sizes for the round jet model was reduced by two-thirds (three-fourths 
for planar jet) to investigate grid dependence of the solution. Since the computational ability was 
limited, overall grid length was reduced (by half) to keep the total number of cells less than 
300,000. All other aspects of the models were kept the same. Figures 67 to 74 in Appendix D 
(Section F) illustrate the results. The plots are shown differently due to the inability of the 
shorter refined model to reach similarity. Therefore, the similarity equations could not be used 


to evaluate the jet characteristics and the predictions had to be displayed as shown. 


The differences between the nominal and refined models were found to be small (<5%), 
therefore grid resolution was deemed satisfactory for both jets. Several additional grid 
configurations were attempted without successful solution convergence. However, future grid 
refinement in the jet shear layer and other high gradient areas may be able to produce more 


significant changes in the solution results. Time constraints limited further research in this area. 
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VII. SUMMARY AND RECOMMENDATIONS 


A critical review of free-jet experimental data was summarized and compared to CFD 
simulations. The CFDRC flow code, using the k-e turbulence model, predicted round 
axisymmetric jets poorly while marginally estimating the two-dimensional planar jet flow. By 
not being able to vary some of the turbulence model constants, the code over predicted the 
velocity spreading and centerline decay of the round jet by approximately 25 percent. Accuracy 
of the planar jet simulations was generally better within an approximately 15 overall percent 


deviation, however larger centerline temperature decay errors existed. 


Sensitivity analysis on both jet simulations indicated that changes in jet exit velocity 
(Ue), turbulence intensity (Tu), and turbulent length scale (L) had little or no effect on the 
solution results. The jet exit temperature (Te) and turbulent Prandtl] Number (Pr,), conversely, 
produced significant changes in the results. The recommended free-jet Prandtl Numbers 
(Chapter III, Section E.2) resulted in correct temperature profiles (spreading rates) when 
compared to the velocity profiles. The change in jet exit temperature also created correct 
centerline decay variations, but failed to produce appropriate spreading rate changes as reported 


in experiments. 


Even though all attempts to increase the accuracy of the simulations failed, the sensitivity 
analysis and insights reported in this study should be extremely useful to follow-on researchers 
utilizing CFDRC’s CFD-ACE+ code. Recommendations for further work include evaluating the 


alternative RNG k-e turbulence model contained within the program. This turbulence model is 


69 


an updated variation of the k-e model that uses a renormalization group approach to 
systematically remove the smallest scales of motion. The RNG k-e turbulence model was tried 
several times without obtaining proper solution convergence, most probably due to grid 
size/refinement limitations (<300,000 cells). | Time constraints prevented further in-depth 


troubleshooting attempts 


More realistic jet exit profiles should also be simulated and compared to experimental 
data. By using Fortran input subroutines within the code, the user can import realistic jet 
velocity and temperature exit profiles. Since jet nozzle boundary layers have a significant effect 
on real jet behavior, their effects and sensitivity on the simulated results should also be 
investigated. In addition, turbulence quantities (k,e) of the initial conditions should be varied to 


investigate their possible effects. 
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APPENDIX A. THE CFDRC PROGRAM 


The CFD Research Corporation (CFDRC) provides a variety of tools for the simulation 
and analysis of fluid flow and associated physics for an assortment of industrial applications. 
CFDRC has specifically developed their software so that the average engineering professional 


can easily manipulate several multi-disciplinary engineering project simulations. 


Typical numerical simulations involve three distinct process steps: 
i 

1) The volume of interest (or solution space) must be divided into discrete control volumes or 
cells. | 

2) One must then define the boundary conditions, initial conditions, and the required equations 
to be solved at each cell. In addition, the numerical technique used to solve the required 
equations must also be defined. 

3) Finally, after the solution has been calculated, the information needed must be extracted from 


the large volume of data generated in the solution process. 


CFDRC provides software modules and packages to address each of these steps for a 
total overall solution. This program, entitled CFD-ACE+ (Version 6.2), uses three separate, yet 
interactive codes to impart an “all in one” seamless commercial code flow solver. CFD-GEOM 
provides interactive geometry modeling and grid generation capabilities. Similar to computer 
aided design (CAD), geometry modeling is the process of creating a computer model of the 


geometry that makes up a problem domain. CFD-GEOM offers a NURBS-based geometry 
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engine with a variety of geometric construction tools. It also has the ability to import and export 
IGES data from most major industry CAD programs. Grid generation is the process of 
discretizing the problem domain with individual cells over which the flow equations are 
integrated. CFD-GEOM provides for the production of two classes of cells: structured and 


unstructured. 


CFD-ACE (U) is an unstructured, polyhedral cell flow solver. It is also integrated with a 
wide variety of physics modules making it the core of a multi-disciplinary analysis environment. 
Inputs are specified for CFD-ACE (U) using CFD-GUI, an advanced graphical user interface that 


allows complete specification of the multi-physics problem. 


CFD-ACE (U) employs a cell-centered control volume solution approach. This approach 
implies that the discrete equations are formulated by evaluating and integrating the fluxes across 
the faces that surround each control volume. In addition, CFD-ACE (U) uses a pressure-based 


methodology in which pressure becomes one of the dependent variables evaluated at each cell. 


The CFD-ACE (U) unstructured flow solver can simulate a wide variety of flow regimes 


and phenomenon as listed below. 
¢ Internal or External Flow 


¢ Laminar or Turbulent Flow: (involves the solution of one or more additional equations) 


* Incompressible or Compressible Flow 


dz 


¢ Heat Transfer 

* Mixing and Reaction: Flows involving multiple gases that mix 
and react can be modeled. 

¢ Steady- state or Transient 


Several others [see CFDRC User Manuals (2000)] 


CFD-ACE (U) also offers a variety of differencing schemes (1*-order upwind to high- 
speed 3-order) to compute the variables within the domain. The type of differencing scheme 
determines how the cell face values are calculated from adjacent cells. These schemes can be 
independently selected for each variable to be solved; however, for this investigation a 1*-order 
upwind scheme was chosen for all variables in the jet models to simplify convergence 


operations. 


CFD-ACE (U) employs an iterative solution technique in which the assembled equations 
for each dependent variable are solved sequentially and repeatedly, with the goal of improving 
the overall solution with each iteration. This solution convergence reduces the variable value 
changes within each control volume until acceptably small values are obtained and the final 
overall solution is determined. The nonlinear coupled nature of the Navier-Stokes (and other 
relevant equations) makes it necessary to restrain or under-relax the iteration-to-iteration changes 
of each variable in order to obtain a stable convergence of the solution procedure. Under- 
relaxation constrains the amount that each variable can change from one iteration to the next. 


The dependent variables (u, v, w, and k) are modified by using an Inertial Factor. However, a 
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linear under-relaxation technique is applied to the auxiliary variables (p, T, 4, 9). The code 


allows the user to change these under-relaxation values within the program to ensure stable 


solution convergence. 


CFD-ACE (U) generates a wide variety of graphical and printed outputs. This output 
includes graphics files with information at each cell for use in CFD-VIEW and the printed 
output, which allows the user to monitor several variables at a fixed location in the problem 
domain. Integrated quantities such as mass flow-rates, heat transfer rates, and pressure forces 
can be written to the text output file. In addition, the change in the solution from iteration to 


iteration, or the residual, can be graphically monitored to assess convergence. 


One of the larger challenges in computational modeling is the management of the large 
volume of data generated with each simulation. These data must be reduced to extract useful 
information, which can be applied to practical problems. To aid in the data reduction process, 


the CFD-ACE+ suite includes a 3-D graphical post-processor called CFD-VIEW. 


CFD-VIEW contains a variety of tools to visualize and extract data from complex 3-D 
data sets. Various types of visual surfaces can be generated such as constant computational 
plane surfaces, cutting plane surfaces, iso-value surfaces, and unstructured surfaces. On each 
surface, contour levels of any of the variables in the data set can be displayed as value-colored 
lines, flooded contour levels, and continuously shaded value-colored surfaces. Vector fields can 


also be displayed on these planes using arrows that indicate magnitude and direction. 
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Other miscellaneous features of CFD-VIEW include a point probe that allows one to 
extract data from any point in the flow field, a line probe used to make X-Y plots, and a 
streamline tracer to create particle traces using any vector field in the model. CFD-VIEW also 
contains an expression calculator that allows one to perform many of the tasks necessary to 
reduce a computational data set into usable information. The calculator has a variety of built-in 
functions that enable cell-wise mathematical operations and derivative mathematical options that 
depend on more than one cell. Results of calculations can be visualized in the same way as any 


other data. 


Overall, the three interactive codes of CFD-ACE+ {CFD-GEOM, CFD-GUI / CFD-ACE 
(U), and CFD-VIEW} work in harmony to step the user through the numerical simulation 


process. The reader is referred to the CFDRC User Manuals (2000) for more comprehensive 


information. 
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APPENDIX B. CFD-ACE (U) STANDARD k-e TURBULENCE MODEL 


The non-linearity of the Navier-Stokes equations, coupled with the complexity of the 
boundary conditions, makes it impossible to obtain analytical solutions for all but a limited 
number of flows of engineering interest. Hence one is forced to resort to approximate or 
numerical methods. Even though a wide variety of numerical techniques can be applied to solve 
the Navier-Stokes equations for laminar flows, Direct Numerical Simulation (DNS) of turbulent 
flows is feasible only at very low Reynolds numbers. Turbulent flows are inherently unsteady 
and they contain a wide range of time and length scales, and resolution of these scales requires 
very short time steps and fine grids. The CPU and memory requirements are too large for most 


present day computers. 


As most engineering applications only require time-mean quantities, the Navier- Stokes 
equations are usually averaged over time or ensemble of statistically equivalent Aaws to yield 
averaged equations. In the averaging process, a flow quantity is decomposed into mean and 
fluctuating parts. Reynolds (or time) averaging and Favre (or density) averaging are two of the 


more popular techniques generally used. 


For turbulent flows, CFD-ACE (U) exclusively utilizes the Favre-Averaged Navier- 
Stokes (FANS) equations to solve for momentum and energy. Although the FANS equations 
contain less information than the full NS equations, they do contain additional unknown 
Reynolds stresses. These correlations between the fluctuating components arise in the averaging 


process, and are additionally modeled to achieve closure of the FANS equations. All the 
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turbulence models available in CFD-ACE (U) employ the generalized Boussinesq eddy viscosity 
concept in which the Reynolds stresses are treated as a linear function of the mean strain rate. 
The Generalized Transport Equation given below indicates the common form for all of the 
FANS equations. 

_%p yat 


Pe ee eases {Source} 


—+5S, (38) 


The symbol @may represent any of the velocity components, enthalpy, or other scalar 
flow variables. In the preceding equation, I'zg is the effective diffusivity and is modeled as: 


i (39) 
ef 
Pr Pr 


Sg is a generalized source term, which represents the mechanisms for the generation and 
destruction of ¢. In addition, any terms that cannot be conveniently expressed as convection or 
diffusion (e. g., the pressure gradient term in the momentum equations) can, in general, be 


lumped into the source term (S4). 


The standard k-e model, employed by CFD-ACE (U), is based on Launder / Spalding 
(1974). The two-equation model, in a FANS generalized form, governs the transport of turbulent 


kinetic energy (k) and its dissipation rate (€). The square root of k is taken to be the velocity 


scale, while the length scale (2) is modeled as: 


f=— (40) 
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The expression for eddy kinematic viscosity is: 


C,k? 
y, = fhe (41) 
p é 
Therefore, the modeled equations for k and € are: 
Q 7) L, \ ok 
k) ee k)= pP- pe+—|| ut+— | 
5, (ok) (on, )= pP-p a, [- za | (42) 
JS per ue Li, \ de 
—C., ae! ac oy Pacts 
a 5 ae k ox, [ "6, iz] (43) 
with the production (P) defined as: 
ou, Ou, 2 Ou, du, 2, Ou,, 
P=v,|—— —*6, | -=—k 
(3 i, Bax, ! fe 3 Si 


The five constants used in the k-€ model are listed below. All of them are internal 


to the code and can not be changed for axisymmetrical flows as suggested by Launder/Spalding 


C,, = 0.09 
C., =1.44 
C,. = 1.92 
0, =1.0 
oO, =1.3 


(1974). Therefore, without being able to vary C, and C,; across the jet profile, the overall 


prediction accuracy for axisymmetrical jets will be poor. 
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The standard k-e model is a high Reynolds Number model and is not intended to be used 
in the near-wall regions where viscous effects dominate the effects of turbulence. Instead, “wall 
functions” are used in cells adjacent to walls. The & and etransport equations are not 
numerically eee in these cells. Instead, semi-empirical expressions are used to relate k, 
¢, and the friction velocity (@). These expressions are obtained from analysis of the momentum 


and turbulence equations for a flat plate boundary layer, assuming a logarithmic velocity profile. 


k= 
oe (45) 
é= (46) 
Ky 
The friction velocity (2) is defined as: 
a= [te 
p (47) 


where tw is the shear stress at the wall and is obtained by assuming the velocity profile between 


the wall and the next-to-wall grid points obeys the following “law of the wall”: 


For y' <11.5 ue=y* (48) 
+ + 1 + 
For y > 11.5 u Sey ) (49) 
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with the dimensionless distances y* and u* defined as: 


yt (50) 
Le 
ut == (51) 
u 


The wall shear stress (tw) is calculated iteratively from the known values of y and w in 
the first cell, The constants appearing in Equations 46 and 49 are experimentally determined to 
be E =9. 0 and x = 0.4 (von-Karman constant). Because the semi-empirical relations for & and 
¢ in the first cell assume a logarithmic velocity profile, the spats wall functions are strictly 
valid only if the center of the cell nearest the wall is inside the logarithmic boundary layer (30 < 
y’ < 150 suggested range). The wall shear stress, evaluated using Equations 48 and 49, is used to 


calculate the boundary condition for the velocity components parallel to the wall. 


For specifying turbulence characteristics at inlet and outlet boundaries, three quantities 
can be utilized; turbulent kinetic energy (A), dissipation rate (€), and length scale (2). The code 
requires that only two of the three quantities be specified. It is usually more convenient, to 
provide a length scale instead of a dissipation rate value. For example, the length scale is usually 
the inlet diameter or height. In this case, the program will internally calculate the boundary 
value of ¢ from the given values of & and ¢ as shown in Equation 52, where the constant, C,, still 


has the value 0.09. 


$1 


© 0.03L (52) 


This formula assumes that the characteristic length scale for the turbulent eddies is 
approximately 0.3% of the size of the inlet dimension [L = D (or H) in our case]. The eddy- 


viscosity (£4, at the boundary is then internally calculated using the given k and the calculated (or 


given) € as: 


2 
ee (53) 


After applying the Favre averaging procedure, the total enthalpy energy equation 


becomes: 


0 0 _ oo oT | dp. a 
ae wiles [Rar tod te, (54) 


t 


where the total (stagnation) enthalpy is defined as: 


h, =it 24+—(u? +v? +0?) (55) 
p 2 
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Similar to the momentum equation analysis, the viscous stress tensor (4%) for Newtonian 


fluids can be related to the velocity gradients by the following equation. 


ei i tn 2 Olle 
aH Oe, | 3) lax | (56) 
where dj is the Kroneker delta function: = §, = 0 fori#j 
% 1 fori=j 
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**Profile (Spreading) Calculations** 
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1.085273 
1.096818 
1.108364 
4.419909 
1.131455 
1.143 


40.0: 
40.07048 


37.05332 


8.537625 
8.398851 
8.243897 
8.093073 
7.958456 
7.822799 
7.681531 
7.563874 
7.435015 
7.313878 
7.184614 
7.077387 
6.963366 
6.851489 
6.744654 
6.635124 
6.53774 
6.441493 
6.342113 
6.246802 
6.164253 
6.077325 
5.99616 
5.913419 
5.84544 
6.772453 
5.707756 


328.4631 
327.8035 


305.7753 
305.6621 
305.5517 
305.4422 
305.3417 
305.2403 
305.1428 
305.0524 
304.96 
304.8752 
304.782 
304.7012 
304.6189 
304,5395 
304.4669 
304.3889 
304.3195 
304.2482 
304.1803 
304.1132 
304.0847 
303.9919 
303.9301 
303.8712 
303.8139 
303.7627 
303.7317 


1.81818189 
2.272727165 
2.727272835 

3.18181811 


33.1818189 
33.6363622 
34.09090945 
34.54545669 


35.45464724 
36,90909056 
36.3636378 
36.8181841 
37.27272835 
37.72727669 
38.1818189 
38,.63636614 
39.09090945 
39,64543307 
40 
40.46462766 
40.90909449 
41.36362206 
41.81818893 
42.27271664 
42,72728346 
43.18181102 
43.63637796 
44.08030551 
44.64647244 
46 


0411241545 
1.024369144 
1.0541428285 
1.079517841 


4.685110906 
4,762522874 
4,852040243 
4,942463759 
S-O2evEess 
6.113223641 
6.207258813 
§.288268371 
6.379911137 
6.469016574 
6.667416246 
6,.66176385 
§.744308141 
6.83810614 
6.930681465 
6.028481156 
6.118279406 
6.209697038 
6.307002099 
6.403231606 
6.488980903 
6.68179709 
6672002082 
&.764228961 
8,842892922 
6.929414583 
7.007968995 


1.000320102 
41.003945506 
1.026479756 
1.053996227 
4.079000845 


5.194535349 
5.298387524 
5§.403750203 
§.512476572 
&.72486308 
6.833398149 
§.937772148 
6.048387097 
6.163593699 
6.273628721 
6381349443 
6.496052935 
6608657341 
8.716087071 
6836428733 
6.946248293 
7.061814416 
7176618432 
7.293691364 
7.398821121 
7.616218317 
7,633393652 
7.749636028 
7.865963974 
7.972998113 
8.039231449 


** Centerline Decay Calculations** 


. 0254 
UMMARY OUTPU rye 


Standard Error 0.008171861 
‘Observations 23 


1 84447753 8.444775385 
21 0.001402365 6.67793E-05 
22 8446177751 


Intercept 
X Vartable 1 


0.99978577 
Standard Error 0.011148643 
bservations 3 


7613985 
21 0.002610137 0.000124292 
22 _12.76400863 


-3.043807526 0.030927481 -98.41757104 1.79884E-29 -3. 
0247048036 0.000771 320.4255852 3.14632E-40 0.2454446§3, 


p' 
X Varlable 4 


APPENDIX D. SENSITIVITY ANALYSIS FIGURES AND TABLES 


JET TURBULENCE INTENSITY 


Table 9. CFD Simulation Predictions with Varying Jet Turbulence Intensity 


C/L Kin. 
Vel. Virt. 
Decay | Orig. 


{Kou} | {Cou} | {Kat 
AXISYMMETRICAL ROUND JET 


{Tu} {Ki} } Ka} | {Cx} 
Target | 0.085- | -1.0to | 0.160- 0.098- 0.185- | 0.0 to 
Values | 0.095 0.172 | 40 | 0110 | 20 | 0.194 | 2.0 


4.0 

0.002 0.200 0.134 0.246 
9.3 
8 


Temp. | Spread 
Spread | Virt. 
Rate Orig. 


{Ci 


Geom. 
Virt. 


{Cru} 


| 0.002_| ; 

[(0.005)on | 0.112 | 9.3 | 0.201 | 100 | 0.134 | i141 | 0.247 | 123 | 
| 0.001__| 
| (0.002)nom_| 
[| 0.01 | 


TW0-DIMENSIONAL PLANAR JET 


Target | 0.096- 0.155- 0.125- 0.250- | -5.0 to 
Values | 0.120 2.0 0.220 0.140 | 0.0 | 0.280 | 00 
0.001 0.091 0.152 | 5.1 | 0.145 8.4 0.207 
u, 9.1 
3 


4.0 
5.1 | 8.4 | 0. 
(0.002),., | 0.092 | 37 | 0.154 | 52 | 0149 | 86 | 0.209 | 9.1 | 
3 8.5 


0.092 | 33 | 0.154 | 
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o Tu=0.01 [K2u=0.203, C2u=9.18] 
Tu=0.005 [K2u=0.201, C2u=10.0] 

A Tu=0.002 [K2u=0.200, C2u=10.2] 
High Target [K2u=0.172, C2u=-1]} 
Low Target [K2u=0.160, C2u=4] 


Ue/Uc 


Figure 27. Round Jet Centerline Vel. Decay with Varying Jet Turb. Intensity 


Oo Tu=0.01 [K2t=0.249, C2t=1 1.4] 
Tu=0.005 [K2t=0.247, C2t=12.3]} 
A Tu=0.002 [K2t=0.246, C2t=12.5] 


(Te-Tinf) / (dTc-Tinf) 


Figure 28. Round Jet Centerline Temp. Decay with Varying Jet Turb. Intensity 


88 


© Tu=0.01 [Kiu= 0.115, Ciu= 8.8] 
Tu=0.005 [K1u= 0.112, Clu= 9.3] 


0.9 3 A Tu=0.002 [Kius.111, C1u=9.4] 
| x High Target [Ktu= 0.095, C1u=-1] 
0.8 a sieeees Low Target [K1u= 0.085, Ciu=4] 
0.7 aR 
"Ty 
i \x'® 
2 0.5 
ne ee 
0 
0 0.05 20.1 015 0.2 0.25 0.3 
Eta=R/(X-CiuD) 
Figure 29. Round Jet Velocity Profile with Varying Jet Turbulence Intensity 
, © Tu=0.01 [Kit= 0.137, C1u= 10.4] 
Tu=0.005 [K1t= 0.134, Ctu= 11.1] 
0.9 A Tu=0.002 [K1t=.134, C1t=11.3] 


(T-Tinf) / (Te-Tinf) 


Figure 30. 


x High Target [K1t= 0.110, C1t=0] 
Low Target [K1t= 0.098, Cit=2] 


0.05 Of O15 02 025 03 
Eta=R/(X-C1tD) 


Round Jet Temp. Profile with Varying Jet Turbulence Intensity 
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Oo Tu=0.001 [K2u=0.152, C2u=5.1] 
Tu=0.002 [K2u=0.154, C2u=5.2] 
4  Tu=0.01 [K2u=0.154, C2u=4.8] 
—--—--High Target [K2u=0.220, C2u=-2] 
Low Target [K2u=0.155, C2u=4] 


(Ue/Uc)*2 


Figure 31. Planar Jet Centerline Vel. Decay with Varying Jet Turb. Intensity 


oO Tu=0.001 [K2t=0.207, C2u=9.1] 
Tu=0.002 [K2t=0.209, C2t=9.1] 
A Tu=0.01 [K2t=0.208, C2u=8.5] 
---.-High Target [K2t=0.280, C2t=-5] 
Low Target [K2t=0.250, C2t=0] 


(Te-Tinf)*2 / (Te-Tint)*2 


Figure 32. Planar Jet Centerline Temp. Decay with Varying Jet Turb. Intensity 
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Tu=0.001 [Kiu=.091, Ciu=3.6] 
Tu=0.002 [K1u=.092, Ciu=3.7] 
4 Tu=0,01 [Kiu=.092, C1u=3.3] 

; x High Target [Kiu= 0.120, C1u=-6] 
0.8 ; ae ae Cee Low Target [K1u= 0.096, C1u=2] 


U/Uc 
eo 
a 


0 O05 O1 O15 O2 025 O03 0.35 
Eta=Y/(X-C1uH) 


Figure 33. Planar Jet Velocity Profile with Varying Jet Turbulence Intensity 


© Tu=0.001 [Kit=.145, Cit=8.4] 

Tu=0.002 [K1t=.149, C1t=8.6] 
A  Tu=0.01 [K1t=.149, C11=8.3] 
x High Target (Kit= 0.140, C1t=-5] 

seceeee Low Target [K1t= 0.125, C1t=0] 


(T-Tinf) / (Te-Tinf) 
-] 
oi 


it) 0.05 0.1 015 O02 025 O03 0.35 
Eta=Y/(X-CitH) 


Figure 34. Planar Jet Temp. Profile with Varying Jet Turbulence Intensity 
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TURBULENT PRANDTL NUMBER 


Table 10. CFD Simulation Predictions with Varying Prandtl Number 


AXISYMMETRICAL ROUND JET 


| CAXISYMMETRICALROUND JET 
Lae 0. 120 4, 0 0. 220 4. 0 0. 140 2. 0 0. 280 2: 0 
Oe foe 3 eae Pee eae 


TW0-DIMENSIONAL PLANAR JET 


[TWO-DIMENSIONAL PLANARJET 
a Ca 0.120 | 120 2. 0 0. 220 4, 0 0. 140 0. 0 0. 280 0. 0 

| 0.091 | 

[O.5)aon | 0.092 |~3.7_| 0.154 | 52 | 0149 | 86 | 0209 | 91_| 


| 07 | 0.0901 | 36 | 0.152 {| 5.1 | 01s | 7.2 [0177 [81 | 
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°o 


A Prt=0.9 [K2u=0.199, C2u=10.0] 
--~--High Target [K2u=0.172, C2u=-1] 
Low Target [K2u=0.160, C2u=4] 


Ue/Uc 


Figure 35. Round Jet Centerline Vel. Decay with Varying Prandtl Number 


© Prt=0.5 [K2t=0.286, C2t=12.7] 
Prt=0.7 [K2t=0.247, C2t=12.3] 
A = Prt=0.9 [K2t=-0.214, C2t=11.5] 


| —.-- High Target [K2t=0.194, C2t=0] 
Pccceeee Low Target [K2t=0.185, C2t=2)} 


(Te-Tinf) / (Te-Tinf ) 
ol 


Figure 36. Round Jet Centerline Temp. Decay with Varying Prandtl Number 
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U/Uc 


Prt=0.5 [Kiu=0.112, C1u=9.3)} 
Prtz0.7 [Kiu= 0.112, Ciu= 9.3] 

A Prt=0.9 [Kiu=0.111, C1u=9.2] 

x High Target [Kiu= 0.095, C1u=-1] 
Low Target [Kiu= 0.085, Ctu=4} 


0 0.05 =O. 0.15 O02 025 03 
Eta=R/(X-C1uD) 
Figure 37. Round Jet Velocity Profile with Varying Prandtl Number 


(T-Tinf) / (Te-Tinf ) 


oO Prt=0.5 [K1t=0.159, Cit=11.4] 
Prt=0.7 [K1t= 0.134, Ciu= 11.1] 
A Prt=0.9 [K1t=0.112, C1t=9.4] 

x High Target [Ktt= 0.110, C1t=0] 
seteeee Low Target [K1t= 0.098, C1t=2] 


0 0.05 0.1 0.15 0.2 0.25 0.3 
Eta=R/(X-C1tD) 


Figure 38. Round Jet Temp. Profile with Varying Prandtl Number 
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© Pr=0.3 [K2u=0.152, C2u=5.1] 
Pr=0.5 [K2u=0.154, C2u=5.2] 

A Pr=0.7 [K2u=0.152, C2u=5.1] 
—-~—--+High Target [K2u=0.220, C2u=-2] 


(Ue/Uc)42 
a 


Figure 39. Planar Jet Centerline Vel. Decay with Varying Prandt! Number 


© Pr=0.3 [K2t=0.240, C2u=8.4] 
Pr=0.5 [K2t=0.209, C2t=9.1] 

A Pr=0,7 [K2t=0.176, C2u=8.1] 
—-—--High Target [K2t=0.280, C2t=-5]} 
a cee Low Target [K2t=0.250, C2t=0] 


(Té-Tinf)*2 / (Tc-Tinf)*2 
on 


Figure 40. Planar Jet Centerline Temp. Decay with Varying Prandtl Number 
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oO Prt=0.3 [K1u=.091, C1u=3.6] 
Prt=0.5 [K1u=.092, C1u=3.7] 
A Prt=0.7 [K1u=.091, Ciu=3.6] 
x High Target [Kiu= 0.120, Ciu=-6] 
veneeee Low Target [Kiu= 0.096, C1u=2] 


Figure 41. 


(T-Tinf) / (Te-Tint ) 


0.1 0.15 O02 0.25 O03 0.35 
Eta=Y/(X-C1uH) 


Planar Jet Velocity Profile with Varying Prandtl Number 


© Prt=0.3 [K1t=.173, C1t=5.6} 
Pit=0.5 [K1t=.149, C1t=8.6] 
A Prt=0.7 [K1t=.117, C1t=7.2] 
x High Target [K1t= 0.140, C1t=-5] 
sees Low Target [K1t= 0.125, C1t=0] 


Figure 42. 


0.1 0.15 0.2 0.25 0.3 0.35 
Eta=Y/(X-C1th) 


Planar Jet Temp. Profile with Varying Prandtl Number 
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C. 


JET LENGTH SCALE 


Table 11. CFD Simulation Predictions with Varying Jet Length Scale 


C, u Kou C, ut Ky t 
AXISYMMETRICAL ROUND JET 


L tu It Ky Cx 
Values 0.095 4.0 0.172 4.0 0.110 2.0 0.194 2.0 
0.199 | 10.0 | 0.133 0.245 
(0127) om | 0.112 | 9.3 | 0.201 | 100 | 0.134 | 111 | 0.247 | 123 | 
0.001 

eee Pale 
| O01 
| (0065) nom_| 
| 0.001 | 


TW0-DIMENSIONAL PLANAR JET 


Values 0.120 2.0 0.220 4.0 0.140 0.0 0.280 0. 
o1 | 0.091 | 3.6 [0152 | 51 [0146 | 84 | 0207 | 9.0 
0.092 | 3.7 {0.154 [52 | o149 | 86 {| 0.209 | 91 | 
3.6 | 0. 5.1 | 84 | 0207 | 9.1 | 


| o.091 | 3.6 | 0.152 | 5.1 | 0.145 | 
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Oo = Le=0.001 [K1u=0.110, C1u=9.2 
Le=0.0127 [Kiu=0.112, C1u=9.3] 
Le=0.1 [K1u=0.111, Ciu=9.1] 

xX High Target [K1u=0,095, C1u=-1] 
sence Low Target [K1u=0.085, C1u=4] 


U/Uc 


0 0.05 60.1 0.15 02 025 03 
Eta=R/(X-C1uD) 


Figure 43. Round Jet Centerline Vel. Decay with Varying Jet Length Scale 


© Le=0.001 [Kit=0.1 
Le=0.0127 [ 


(T-Tinf) / (Te-Tinf) 


0 0.05 O01 O15 O02 025 03 
Eta=R/(X-C1tD) 


Figure 44. Round Jet Centerline Temp. Decay with Varying Jet Length Scale 
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Le=0.001 [K2u=0.199, C2u=10.1] 
Le=0.0127 [K2u=0.201, C2u=10.0] 
A Le=0.1 [K2u=0.199, C2u=10.0] 
-—- - High Mus acre C2u=-1]) 
ceeeees Low Target [K2u=0.160, C2u=4] 


Ue/Uc 
oo 


Figure 45. Round Jet Velocity Profile with Varying Jet Length Scale 


o Le=0.001 [K2t=0.244, C2t-12 
Le=0.0127 [K2t=0.247, C2t=1 

A Le=0.1 [K2t=0.245, C2t=12.3] 
—-~—--High Target [K2t=0.194, C2t=0) 
seeeeee Low Target [K2t=0.185, C2t=2] 


(Te-Tinf) / (Te-Tint) 
a 


Figure 46. Round Jet Temp. Profile with Varying Jet Length Scale 
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(Ue/Uc)*2 


oO Le=0.001 [K2u=0.152, C2u=5.1] 
Le=0.0065 (K2u=0.154, C2u=5.2] 

A Le=0.1 [K2u=0.152, C2u=5.1] 
—----High Target [K2u=0.220, C2u=-2] 
Low Target [K2u=0.155, C2u=4] 


Figure 47. Planar Jet Centerline Vel. Decay with Varying Jet Length Scale 


(Te-Tinf)42 / (Tc-Tinf)*2 


Figure 48. 


o = Le=0.001 [K2t=0.207, C2u=9.1] 
Le=0.0065 [K2t=0.209, C2t=9.1] 
A Le=0.1 [K2t=0.207, C2u=9.0] 
—-—--High Target [K2t=0.280, C2t=-5] 
Low Target [K2t=0.250, C2t=0] 


Planar Jet Centerline Temp. Decay with Varying Jet Length Scale 
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o  Le=0.001 [Kiu=.091, C1 uae) 
Le=0.0065 [K1u=.092, C1u=3.7] 
A Le=0.1[Ktu=.091, C1u=3.6] 
x High Target [Kiu= 0.120, Clu=-6] 
wees Low Target [K1u= 0.096, Ciu=2] 


U/Uc 
fo] 
oi 


0 005 O1 O15 O2 O25 03 0.35 
Eta=Y/(X-C1uH) 


Figure 49. Planar Jet Velocity Profile with Varying Jet Length Scale 


o  Le=0,001 [Kit=.145, C1t=8.4] 
Le=0.0065 [K1t=.149, C1t=8.6] 

A Le=0.1[K1t=.146, C1t=8.4] 

x High Target [Kit= 0.140, C1t=-5] 
Low Target [K1t= 0.125, C1t=0] 


(T-Tinft) / (Te-Tinf ) 


0 005 O1 O15 O02 025 O03 0.35 
Eta=Y/(X-C1tH) 


Figure 50. Planar Jet Temp. Profile with Varying Jet Length Scale 
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D. 


JET VELOCITY 


Table 12. 


Values 


0.085- 
0.095 


-1.0to / 0.160- | -1.0to } 0.098- | 0.0to | 0.185- | 0.0 to 
4.0 0.172 4.0 0.110 2.0 0.194 2.0 


CFD Simulation Predictions with Varying Jet Velocity 


{Cou} | {Kae} | {Cie} 
AXISYMMETRICAL ROUND JET 


0.112 9.3 0.201 0.134 0.247 
0.113 


Values 0.120 2.0 0.220 4.0 0.140 0.0 0.280 
9.1 
| 86 


TW0-DIMENSIONAL PLANAR JET 


9.4 
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oO Ue=20 [K2u=0.194, C2u=9.5] 
Ue=40 [K2u=0.201, C2u=10.0] 

A Ve=60 [K2u=0.203, C2u=10.2] 
— - -High Target [K2u=0.172, C2u=-1] 
By eoedsts Low Target [K2u=0.160, C2u=4] 


Ue/Uc"” 
a 


Figure 51. Round Jet Centerline Vel. Decay with Varying Jet Velocity 


@  Ue=20 [K2t=0.250, C2t=12.4] 
Ue=40 [K2t=0.247, C2t=12.3] 

A Uex60 [K2t=0.240, C2t=12.2] 
—-—--High Target [K2t=-0.194, C2t=0) 
OF), sectors Low Target [K2t=0.185, C2t=2] 


(Te-Tinf) / (Tc-Tinf) 
oa 


Figure 52. Round Jet Centerline Temp. Decay with Varying Jet Velocity 
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°o 


Ue=20 [K1u=0.111, C1u=9.1) 
Ue=40 [K1u=0.112, C1u=9.3] 
A  Ue=60 [Kiu=0.112, C1u=9.4] 
x High Target [K1u=0.095, C1u=-1]} 
teens Low Target (K1u=0.085, C1u=4] 


U/Uc 


0 0.05 0.1 0.15 0.2 0.25 0.3 
Eta=R/(X-CiuD) 


Figure 53. Round Jet Velocity Profile with Varying Jet Velocity 


Oo Ve=20 [K1t=0.130, C1t=10.8] 
Ue=40 [K1t=0.134, Ciu=11.1] 

A Ue=60 [K1t=0.136, Cit=11.1} 

x High Target [K1t=0.110, C1t=0) 
Low Target [K1t=0.098, C1t=2] 


(T-Tinf) / (Tc-Tinf) 


0 0.05 0.1 0.15 0.2 0.25 0.3 
Eta=R/(X-C1tD) 


Figure 54. Round Jet Temp. Profile with Varying Jet Velocity 
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Oo Ue=15 [K2u=0.153, C2u=5.1] 
Ue=35 [K2u=0.154, C2u=5.2] 

A Ue=55 [K2u=0.152, C2u=5.1] 
~-—--~High Target [K2u=0.220, C2u=-2] 


ceeneee Low Target [K2u=0.155, C2u=4] ra 

ie ae a ae 
Pt TT et a 
fat aL ae 
oe ee ee a 


(Ue/Uc)42 
oi 


Figure 55. Planar Jet Centerline Vel. Decay with Varying Jet Velocity 


o Ue=15 [K2t=0.212, C2u=9.1] 
Ue=35 [K2t=0.209, C2t=9.1] 
A Ue=55 [K2t=0.198, C2u=8.8] ||, 
~.—--High Target [K2t=0.280, C2t=-5] 
84....... Low Target [K2t-0.250, C2t=0] 


(Te-Tinf)*2 / (Te-Tinf)42 
a 


Figure 56. Planar Jet Centerline Temp. Decay with Varying Jet Velocity 
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“U/Uc 


(T-Tinf) / (Tc-Tinf ) 


Figure 57. 


Figure 58. 


Ue=15 [K1u=.091, C1u=3.5] 
Ue=35 [K1u=.092, C1u=3.7] 

4 Ue=55 [Kiu=.091, C1u=3.6] 

x High Target [K1u= 0.120, C1u=-6] 


0.1 015 O2 0.25 O03 0.35 
Eta=Y/(X-C1uH) 


Planar Jet Velocity Profile with Varying Jet Velocity 


Oo Ue=15 [K1t=.146, C1t=8.3] 
Ue=35 [K1t=.149, C1t=8.6] 

A Ue=55 [K1t=.146, C1t=8.4] 

x High Target [K1t= 0.140, Cit=-5] 
Low Target [K1t= 0.125, C1t=0] 


0.2 
Eta=Y/(X-C1tH) 


Planar Jet Temp. Profile with Varying Jet Velocity 
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E. JET TEMPERATURE 


Table 13. CFD Simulation Predictions with Jet Temperature 


Geom. in. | Temp. | Spread 


Virt. Spread | Virt. 
Orig. 


{Cyu} 
1.0 to 


lu 
AXISYMMETRICAL ROUND JET 


Values 0.095 4.0 0.172 4.0 0.110 2.0 0.194 2.0 
F500 | 0112 | 92 | 0.244 | 104 | 0.133 | 108 | 0.309 | 126 | 


TW0-DIMENSIONAL PLANAR JET 


Values 0.120 2.0 0.220 4.0 0.140 0.0 0.280 0.0 
F G30)oom | 0.092 | 3.7 | 0.154 | 5.2 | 0149 {| 86 | 0.209 | 9.1 
soo | 0.093 | 42 | 0.148 [ 384 | 0.133 | 10. | 0.301 | 10.5_| 


~— 
enya? 
_— 
Le 
[7 
a] 
OQ 
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Te=330 (K2u=0.201, C2u=10.0] 
A Te=500 [K2u=0.244, C2u=10.4] 
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Figure 59. Round Jet Centerline Vel. Decay with Varying Temperature 
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Figure 60. Round Jet Centerline Temp. Decay with Varying Temperature 
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Figure 61. Round Jet Velocity Profile with Varying Temperature 
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Figure 62. Round Jet Temp. Profile with Varying Temperature 
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Figure 63. Planar Jet Centerline Vel. Decay with Varying Temperature 
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Figure 64. Planar Jet Centerline Temp. Decay with Varying Temperature 
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Figure 65. 
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Figure 66. 
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Planar Jet Velocity Profile with Varying Temperature 
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Planar Jet Temp. Profile with Varying Temperature 
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F. GRID REFINEMENT 
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Figure 67. Round Jet Centerline Vel. Decay with Grid Refinement 
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Figure 68. Round Jet Centerline Temp. Decay with Grid Refinement 
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Figure 69. Round Jet Velocity Profile with Grid Refinement 
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Figure 70. Round Jet Temp. Profile with Grid Refinement 
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Figure 71. 
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Figure 72. 
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Planar Jet Velocity Profile with Grid Refinement 
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